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We present numerical results for 2-d CP(N-l) models at 8 = and n obtained in the D-theory formulation. 
In this formulation we construct an efficient cluster algorithm and we show numerical evidence for a first order 
transition for CP(N — 1>2) models at 9 = tt. By a finite size scaling analysis, we also discuss the equivalence in 
the continuum limit of the D-theory formulation of the 2-d CP(N — 1) models and the usual lattice definition. 



1. Introduction 

CP(N — 1) models £Q are defined on the man- 
ifold S 2N ~ 1 /S 1 which is the coset space of the 
breaking of SU(N) down to U(N — 1). The field 
P{x) is parametrized by Hermitean NxN projec- 
tor matrices, i.e. P(x) 2 = P(x) and Tr P(x) = 1. 
The action is given by 

= ^J dx WW^Wl. (!) 

and it is invariant under the global f2 £ SU(N) 
transformation P(x)' = QP(x)Q + . The 2-d 
CP(N— 1) models are interesting toy-models of 4- 
d Yang-Mills theories. In fact, these two quantum 
field theories share several important features, in- 
cluding asymptotic freedom, the dynamical gen- 
eration of a mass-gap and an instanton topolog- 
ical charge. In particular, this latter property 
leads to a non-trivial vacuum angle 9 and allows 
to consider CP(N — 1) models in order to inves- 
tigate ^-vacuum effects. It has been conjectured 
that CP(N — 1) models have a phase transition 
at 6 — tt . For CP {I) = 0(3) this phase transi- 
tion is known to be second order with a vanish- 
ing mass-gap |2I3I4| . For N > 2 the transition 
is conjectured to be first order [5I6| . Studying 
these nonperturbative problems is highly nontriv- 
ial and one can eventually address these questions 
only numerically. However, the numerical sim- 
ulation of CP(N-l) models is not straightfor- 
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ward, even at 6 = 0. In fact, in contrast to O(N) 
models where the Wolff cluster algorithm [7] can 
be used, no efficient cluster algorithm is avail- 
able | 8I9 |. Still a rather efficient multi-grid algo- 
rithm was developed in |10| . At 9 = tt the numer- 
ical investigations become even more difficult due 
to a very severe sign problem. In fact, in this case, 
the topological sectors with even and odd charges 
have opposite Boltzmann weights and their con- 
tributions almost completely cancel in the parti- 
tion function. This makes it exponentially hard to 
simulate the large lattices necessary to reach reli- 
able conclusions on the phase structure. For this 
reason, efficient simulations at 9 = tt have been 
performed only in the CP(1) = 0(3) case by a 
Wolff- type meron-cluster algorithm J]. Here we 
present results obtained with a method allowing 
to perform an unbiased and very accurate numer- 
ical study of any 2-d CP(N - 1) model at 9 = 
and 7T Our method is based on the D-theory 
formulation of a field theory ^21- This is a dif- 
ferent formulation than Wilson's and consists in 
obtaining the continuous fields from the dimen- 
sional reduction of discrete variables. It provides 
an alternative lattice regularization of field theory 
but with the same continuum limit. More details 
on the results presented here can be found in |11| . 

2. CP(N-l) models: D-theory formulation 

In this section we discuss 2-d CP(N-l) models 
in the D-theory formulation. The discrete vari- 
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ables we consider are SU (-/V) quantum spins — 
\X%, generators of the SU(N) group [T" , Ty] = 
idxyfabcTx- The spins are located on the sites x 
of a I x L' periodic square lattice with L 3> £' • 
Hence we have the geometry of a spin ladder 
consisting of n = L'/a transversely coupled spin 
chains of length L. We consider n.n. couplings: 




Figure 1. Spin ladder geometry: open circles form 
sublattice A, filled circles belong to sublattice B. 
antiferromagnetic along the chains and ferromag- 
netic between different chains. The lattice nat- 
urally decomposes into two sublattices A and B. 
The spins in A transform in the fundamental rep- 
resentation {N}, while those belonging to B are 
in the conjugate one {N} and are described by 
—T£*. The quantum spin ladder Hamiltonian is 



xEA 



(2) 



xEB 



where J > 0, and 1 and 2 are unit-vectors in the 
x- and y-directions, respectively. The system has 
a global SU(N) symmetry, i.e. [H, T a ] = 0, with 
the total spin T a = J2 xeA T x ~ Exeb T **> satis- 
fying the SU(N) algebra [T a , T b ] = if abc T c . The 
quantum partition function at the temperature 
T = 1/(3 is then given by Z = Tr exp(-/3if). 

Using the coherent state technique of ^2) one 
finds that, at T = and for L, L' — > oo, the 
SU(N) symmetry breaks down to Z7 (JV— 1). Thus 
there are massless Goldstone bosons described by 
fields PeSU(N)/U(N-l) = CP(N-l). When 
the y-direction is compactified to a finite extent 
L' , the Coleman-Hohenberg-Mermin- Wagner the- 
orem forbids massless excitations. Thus the Gold- 
stone bosons pick up a nonperturbatively gener- 
ated mass-gap m = l/£ and the correlation length 
£ is now finite. For sufficiently many transversely 
coupled chains, the correlation length becomes 
exponentially large, £ 3> L', and the system both 
undergoes dimensional reduction and approaches 



the continuum limit. Using chiral perturbation 
theory, one finds that the lowest-order terms in 
the 2-d Euclidean effective action are given by 

S[P] = J dxdt{-^ Tv^Pd^P] - i6Q[P}}. (3) 

where Q[P] is the topological charge 

Q[P] = ^~J dxdt ^TriPd^PduP], (4) 

and 9 — mi is the vacuum angle. Hence, 9 = 
for even n and 9 = n for odd n. 

3. The numerical study 

One advantage of D-theory is that it is a for- 
mulation of a quantum field theory in terms of 
simple discrete degrees of freedom instead of the 
usual continuum classical fields. In particular, 
the partition function Z of the SU (N) quantum 
spin ladder can be written using a basis of dis- 
crete SU(N) spin states: q £ {u,d, s, ...} on sub- 
lattice A and q £ {u, d, s, ...} on sublattice B. 
These can be simulated with the loop-cluster al- 
gorithm |14I15| . The efficient cluster algorithm 
we have used in this study is a generalization to 
SU(N) of the SU(2) case Q2]. 

The Hamiltonian (J2J) is the sum of 4 terms that 
we denote, respectively, with Hi, . . . , H^: two of 
them (7?2, H4) describe the ferromagnetic inter- 
action while the other two (Hi, H^) the antifer- 
romagnetic one. Using the Trotter decomposition 
we write Z as follows 



Z = lim Trie""" "<• 



(5) 



where we have introduced AM slices in the /In- 
direction with spacing e = (3J/M. Inserting 
complete sets of states between the exponential 
factors, one obtains that the only non- vanishing 
transfer matrix elements are given by 

(6) 



(q\q l \e~ 


~ eH2 |gW) = e £ ^ 




( 




e- tH2 \q l ,q^ 1 ) = (e £ ^ 


+ e 




(q\q^ 
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N + l 

£ —)/2 


(q\q l \e 


- eHl \q l ,q l ) = (e* + iV - 
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-zk /N 
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The transfer matrix elements for H 3 and H4 can 
be obtained by substituting Hi — > H 3 , H 2 — > H4 
and q^q. With these matrix elements at hand, 
the construction of the quantum spin cluster and 
its updating proceeds like in the SU(2) case. One 
can also avoid the slicing in the /3-direction oper- 
ating directly in continuous Euclidean time |16| . 

We have used this cluster algorithm to investi- 
gate if CP(N - 1) models with N > 2 have a first 
order phase transition at 9 = tt where the charge 
conjugation symmetry C is spontaneously bro- 
ken. An order parameter for such a phase transi- 
tion would be the topological charge Q[P] which 
indeed is C-odd: Q[ c P] = Q[P*} = -Q[P]. How- 
ever Q is defined only in the framework of the tar- 
get continuum theory. In the discrete spin system 
we define as order parameter Q[q, q] the number 
of spin flips in a configuration. Q[q, q] gets a con- 
tribution + 1 if a pair of n.n. spins along the re- 
direction, q x q x+ i, flips to another state, q'xrf.p 
at some moment in time. A spin flip from q x q x+ i 
to q* x q'^ instead contributes — 1. In the quan- 
tum spin ladder, charge conjugation corresponds 
to replacing each spin state q x by q x+ \ and, in- 
deed, Q[q,q] changes sign under this operation 
while the action remains invariant. 

We have simulated SU(N) quantum spin lad- 
ders for N = 3,4, and 5, with n = L'/a e 
{2, 3,..., 7}. Our simulations confirm the exis- 
tence of a first order phase transition with spon- 
taneous C-breaking at 9 — ir for all N > 2. As ex- 
pected, there is no phase transition at 9 = 0. Fig- 
ure shows Monte Carlo time histories of Q[q,q] 
for ST/ (4) spin ladders with n = 3 and 4. For 
n = 3 one observes two coexisting phases with 
spontaneous C-breaking, while for n = 4 there is 
only one C-symmetric phase. Similar results have 
been obtained in the other cases studied. 

Finally, we have verified explicitly that the 
CP(N — 1) models in Wilson's and in the D- 
theory formulation have the same continuum 
limit. In order to do this we have compared a 
physical quantity - namely the universal finite- 
size scaling function F(£(L)/L) = f(2L)/£(L) - 
measured in both frameworks. In we show 
the excellent agreement between the results ob- 
tained in these two formulations. 




Figure 2. Monte Carlo time histories of Q[q,q] 
for the CP (3) model at 9 — and ir. 
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